arXiv:1506.00964vl [cond-mat.mes-hall] 2Jun2015 


1 


General Geometric Fluctuation Modeling for 
Device Variability Analysis 

Bo Fu, Seonghoon Jin, Woosung Choi, Keun-Ho Lee, Young-Kwan Park 


Abstract 

The authors propose a new modeling approach based on the impedance field method (IFM) to analyze the general 
geometric variations in device simulations. Compared with the direct modeling of multiple variational devices, the 
proposed geometric variation (GV) model shows a better efficiency thanks to its IFM based nature. Compared with 
the existing random geometric fluctuation (RGF) model where the noise sources are limited to the interfaces, the 
present GV model provides better accuracy and wider application areas as it transforms the geometric variation into 
global mesh deformation and computes the noise sources induced by the geometric variation in the whole simulation 
domain. GV model also provides great insights into the device by providing the effective noise sources, equation-wise 
contributions, and sensitivity maps that are useful for device characterization and optimization. 


I. Introduction 

Device variability presents a growing challenge as we push forward Moore’s Law. Various physical variability 
sources such as random dopant fluctuations (RDF), random workfunction fluctuations (RWF), and random geometric 
fluctuations (RGF) are becoming important as devices get smaller 0, Q, 0, f4j, 0. Among these variability 
sources, geometric variation becomes more significant as the uncertainty in the fabrication process such as lithogra¬ 
phy, etching, and deposition cannot be scaled with the advancing technology nodes. For example, line edge roughness 
(LER) caused by chemical inhomogeneity in the lithography and etching process is in the order of several nanometers 
and is not expected to decrease as the device size shrinks El. Geometric variation induced device fluctuations in 
threshold voltage, drive current, subthreshold swing, and drain-induced barrier lowering (DIBL) exhibit significant 
challenges for future VLSI circuit technology development. 

In terms of the state of the art geometric variation analysis, people either apply direct modeling with sufficient 
device samples 0, 0 or RGF model 0. The former is reliable, however it greatly suffers in terms of efficiency. 
Thanks to the impedance field method (IFM) in, the latter can provide efficient solutions if variations are small 
such that they are in the linear response range. However, since the RGF model constrains the noise source to be 
locally at the interface, it fails to consider nonlocal effects in geometric fluctuations. It is worth noting that the 
recent improvement of RGF model was made to capture the semiconductor body thickness variations. However, 
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the accuracy is not guaranteed, and it relies on the calibration of parameters 0. The purpose of this study is to 
propose a rigorous IFM based model named GV (geometric variations), calculating general geometry variations in 
an efficient and accurate manner. 

The paper is organized as follows: RGF model’s limitations will be explored first regarding two geometric variation 
types: interface displacement and body deformation. Then we explain in detail the formulation of our GV method. 
We will then validate our model’s effectiveness by comparing with numerical experiments on different geometric 
variation cases in the two dimensional device simulations. For a three dimensional application, we demonstrate the 
GV model’s usage on a FinFET structure in which we have considered channel length, body thickness, and fin 
height variations. During the process of the FinFET GV simulations, new physical insights have been obtained such 
as the effective noise sources, equation-wise contributions, and sensitivity maps. Finally, we make a summary and 
conclusion. 


II. RGF Limitations 

We have tested Sentaurus™ Device RGF model 0 for interface displacement variability calculations. However, 
inconsistent results are observed in comparison with the designed numerical experiments. By checking the calculated 
effective Green’s functions for RGF and examining the noise source described in 0, we conclude that the RGF 
model is incapable of modeling general geometric variations. 

A. Interface Displacement 

RGF model accounts for the interface displacement and its low frequency noise source is given as ED: 

K(fi,fi) = n(r~i) ■ n(f* 2 )exp(- — ) 

[^iso (r!) + \a(fi) • n(r"i)|] 

[fliso(^2) + 1^2) • n(r 2 )\], (1) 

where f\ and rj are positions on the interface, n(f) is the local interface normal in position f, a{ S0 the isotropic 
correlation amplitude, a(r) the vectorial correlation amplitude, and A is the correlation length. 

To test the effectiveness of the model, a simple MOSFET structure as shown in Fig. |TJb) is designed. The device 
is composed of only silicon and oxide, with source/drain donor concentration of 3.0 x 10 20 cm -3 and acceptor 
concentration of 3.0 x 10 16 cm -3 in all silicon region. We applied the RGF model on the “silicon/oxide” interface 
that is normal to the Y axis. By choosing a sufficiently large correlation length A = 1 m, the exponential term in 
0 vanishes. Physically it is equivalent to shifting the whole interface simultaneously. Further by setting a(f) = (0, 
0.1 nm, 0) and a{ S0 = 0, the displacement is fixed to be 0.1 nm along the Y-axis direction. Accordingly, we can 
create a variational device as shown in the right of Fig. |TJc) for the numerical experiment, in which the interface is 
moved from silicon to oxide by 0.1 nm. Numerical experimental results are obtained by calculating the difference 
of drain currents between the interface-shifted device and the nominal device. 
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Fig. 1. (a) Current Green’s function for geometric variation, which only has non-zero values at the interface, (b) The nominal device, whose 

sizes are as shown in the X-Y coordinates, (c) An illustration of interface displacement. Whole silicon/oxide interface is shifted from y = 4 
nm to y = 4.1 nm. Only the white circle area in (b) is zoomed for the nominal device (left) and variational device (right). 


Noise-like IFM for RGF calculates the autocorrelation noise current spectral densities (SI). Giving our noise 
source K (r~i, rj) = (0.1 nm) 2 along the Y-axis, we can compare the square root of SI with the absolute values 
of numerical experimental results regarding drain current variations, as shown in Fig. [2] Philips unified mobility 
model iflOl is used in Fig.|2ja) and enhanced Lombardi model ifTTIl . lfl2l is added in Fig.^b). Obvious discrepancy 
exists in Fig. |2jb) because the Lombardi model accounts for mobility degradation near the semiconductor/insulator 
interface, while RGF only captures the noise source locally at the interface as pointed out by Q and demonstrated 
by the effective Green’s function in Fig. [lja) which exhibits non-zero values only at the interface. Specifically, in 
the presence of Lombardi model, the low field mobility is given by ifTO . (T2]: 

— =-b D[ - H-], (2) 

V V b Mac(^±) Msr(^J_) 

where fi^ is the bulk mobility, /i ac and /i sr are the surface acoustic phonon mobility and surface roughness mobility 
both depending on the transverse electric field normal to the interface F_\_, the damping factor is defined as 0: 

D = exp(— -^— ), (3) 

Lrit 

where l CT i t is a fitting parameter with a default value of 10 nm, and l is the distance from the interface. Interface 
displacement changes the l and F± values of the nodes that are not at the interface, and hence their mobilities. 
This effect is not considered in the RGF model. 

Similarly, RGF model fails when thin-layer mobility model lfl3l is applied. For high-field saturation model fl4l . 
m, the driving force of the node not at the interface can be also changed as it is proportional to the gradient of 
scalar variables. By considering just the noise source locally at the interface, RGF is also prone to give inaccurate 
results. 
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Gate Voltage (V) Gate Voltage (V) 

Fig. 2. Results comparison of Sentaurus™ Device RGF model with designed numerical experiments as described in Fig. [7] Philips unified 
mobility model is applied in (a). In addition to the Philips unified mobility, enhanced Lombardi model is added in (b). 


B. Body Deformation 

Besides interface displacement, more general geometric variation is body deformation. For the MOSFET in 
Fig. |TJb), body deformation may include but is not limited to silicon width, oxide thickness, and channel length 
variations. The recent improvement of RGF model added three options of calculating the thickness variation l8l . 
Among the three options, one is used to control the calculation performed on one side or both sides of the interface, 
and the other two are fitting parameters. The basic assumption of the RGF model remains the same, which only 
relies on the local interface information as noise source to predict terminal variations. Thickness variation could 
be introduced by properly setting interface displacement in the latest RGF model. However, RGF model requires 
users to calibrate the two fitting parameters, which demonstrates that the physics is not fully captured. It becomes 
worse if we want to apply the RGF model for length variation because there is no clear and reasonable interface 
to pick as the origin of the noise source. 


III. GV Methodology 

A. GV Formulation 

For continuum device simulations, we solve the following three coupled equations: 

• Poisson’s equation 

J 7 * := -Ve • V<f> - p = 0, (4) 

where e is the permittivity, <f> is the electrostatic potential and p is the net charge. 

• Electron and hole continuity equations 

Tj := V-F + # = 0, (5) 

where F is the carrier (electron or hole) flux, and R is the net recombination rate. 
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Fig. 3. An illustration of finite volume method discretization in 2D. For nodes i and j, Wij is the cross-section length and Lij is the 
side-length. (dashed blue box) is the control volume of node i. Geometric variation moves node j closer to node i. 


Electron and hole density gradient quantum potential equations nsi. For the electron case 

7 h 2 V 2 v^ 


'■= 


- An = 0, 


6m yjn 

where 7 is a fitting parameter, m is the density-of-states effective mass, and A is the quantum potential, 
can be also solved in terms of the quantum potential formulation fI71 . 

A general integral form of each equation a in (<E>, J, A) is given by: 


( 6 ) 

© 


M a := (p F • nd^f - I Sdcr = 0, (7) 

Jr Jn 

where F and S are the flux and source terms respectively, T is the boundary of control volume Q and n is the 
normal vector. 

After obtaining the steady-state solution for the nominal device with mesh M, we can calculate the current 
Green’s Function G l a for terminal i and equation a using the generalized adjoint method m, 03- Since voltage 
and current Green’s functions are mutually transformable using the impedance and admittance matrices, we only 
focus on the current calculation. The calculation of terminal voltage variation follows in a similar way. 

Suppose the geometric variation introduces a new mesh M', whose node coordinates are different from M but 
topologically M' is equivalent to M, i.e. the number of nodes and their connectivity remain unchanged. M' can 
be used to reassemble |4])-([6]) through ([7]) with the already available steady-state scalar solutions (0, n, p , A n , 
A p , etc.) from the nominal device. For equation a , the difference between the RHS on M' and the RHS on M, 
5M a = M a (M f ) — M a {M ), is ultimately the effective noise source due to the geometric variation. The terminal 
current change then can be calculated as follows: 


sr = 

A) 

summing the contribution from each equation. 


( 8 ) 


B. Effective Noise Source Calculation 

Finite volume method (FVM) f20l can be used to solve Q- A plot of FVM discretization is illustrated in Fig. [3| 
in which box coefficient can be defined as Aij = Wij/Lij, the cross-section to side-length ratio, and node volume 
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U is the area of the dashed blue box. Intuitively, geometric variation changes only box coefficient A and node 
volume Vt. Under this assumption, the effective noise source from equation a can be obtained as: 


SM a = 


m °LSA+^6Sl 


dA 


on 


(9) 


where 5A and SQ are the box coefficient and volume changes, respectively. This is similar to the method by A. 
Gnudi, et al CD- Instead of applying Green’s function on each equation as in ([8]), their approach is to convert (5Q, 
SA) to (S(/) 9 Sn , Sp) to continue the variability calculations. 

However, the assumption that only 5A and 5VL change is not always valid. Three additional consequences of 
geometric variations should be considered: 

• Shift of node j as shown in Fig. [3] changes all the distances between j and other nodes. If a semiconduc¬ 
tor/insulator interface consists of node j, a different distance value affects the flux term at the other node. This 
component is critical to correctly simulate geometric dependence of, for example, the enhanced Lombardi 
model. 

• Within the elements where node j belongs to, element gradient coefficients will be different, which alters 
element vectors such as electric field calculated from scalars. Geometric dependence of high field mobility 
usually requires this component. 

• If node j happens to be on the internal boundary, a proper boundary condition should also be taken into 
account, which occurs for example when we want to simulate wave function penetration into insulator using 
([6]) by applying an inhomogeneous Neumann boundary condition fTTl . 

In order to emphasize the importance of all components, we have identified and extracted them during reassem¬ 
bling the RHS of 0 on the variational mesh M'. Detailed comparison of each component’s contribution is shown 
in the application section. 


C. Geometric Variation Mesh Generation 

The present GV model requires that the variational mesh M' and nominal device mesh M share the same 
topology. For this purpose, we have created a strained mesh scheme to only change vertex coordinates. In three 
dimension, for vertex (xq 9 yo, z$) 9 the change of coordinate (5xo,Syo,5zo) is: 

5xq = Ax • fxx(x 0 ) • fxY (?/o) * fxz(zo) 

< Sy 0 = Ay • /yxO 0 ) • /yy(2/o) • ZyzOo) ( 10 ) 

Szo = Az • /zx(^o) * /zy(j/o) • fzz(zo) 

where ( Ax , Ay, Az) are predefined amplitudes of geometric change along the (X, Y, Z) directions. A modulation 
function, for example /xyG/)» is defined on the X-Y plane, describing the y 9 s influence on x. Consequently, the 
vertex (xo, yo, zo) in nominal device mesh M will be moved to (xo + yo + Syo, zo + £zo) in the new mesh 
M'. 
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(a) 



(b) 



Fig. 4. A toy application of strained mesh scheme to mimic geometric variations: (a) a square is transformed to have (b) inward sine functions 
on four sides and (c) outward sine functions on four sides. All shapes have the same node connectivity but varied node coordinates. 


Fig. [^demonstrates a toy application of the strained mesh scheme: every side of a square (Fig.|4ja)) is modified by 
the inward sine function (Fig.[4jb)) or by the outward sine function (Fig.[4jc)). Although shapes change significantly 
from Fig. |4ja) to (b) and (c), their node connectivities remain the same. More practical applications of the strained 
mesh schema will be shown in the application section. 

IV. Applications 

Using our in-house device simulation tool [22] which includes the GV model, we have performed the following 
device simulations. Regarding mobility models, we have employed Philips unified mobility Col, enhanced Lombardi 
surface mobility fTO . lfT2l . and high field saturation mobility lH4l . fl5l . 

The effectiveness of our GV model can be validated by comparing with numerical experiments. Consider that the 
strained meshes created for GV analysis are variational devices that we can run direct simulations on. Drain current 
change A/ D induced by a geometric variation therefore can be obtained by comparing the Ij^-Vq characteristic of 
a variational device on M' and the nominal device on M, i.e. A/d = In(M') — Jd(M), which should agree with 
the result calculated by the GV model. 

For the two dimensional device application, we mainly focus on the model validation by studying both interface 
displacement and body deformation. In the practical three dimensional FinFET example, based on comparison 
with numerical experiments, linear range of GV model will be examined. Contributions of effective noise source 
components besides box coefficients will be compared. We also look at the equation-wise contributions and the 
geometric sensitivity maps. 

A. Two Dimensional DG-FET 

Consider a double gate field effect transistor (DG-FET) as shown in Fig. [5] All dimensions are marked in the 
plot. The current Green’s functions for ([4])-([6]) are drawn in Fig. [6] when both the drain and gate voltages are at 
0.8 V. We will validate the GV model with numerical experiments examining silicon/oxide interface displacements 
and silicon body deformations in the channel region. 

1) Interface Displacement: Depending on the correlation between the two gates, (a), (b), and (c) in Fig. [7] show 
three different cases of silicon/oxide interface displacements. The maximum thickness change of the gate interfacial 
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Fig. 5. Simulated 2D DG-FET composed of silicon, SiC> 2 , FEO 2 , and Nitride. Plotted donor concentration is in the log scale and a uniform 
acceptor concentration of 3.0 x 10 18 cm -3 is used in the simulation. 
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Fig. 6. Current Green’s functions: (a) Poisson’s equation (b) electron continuity equation (c) electron density gradient equation. Vd=Vg=0.8 
V. 


oxide in the middle of the channel is 0.5 nm at each side (equals to 12.5% width variation) which is a large value 
for the illustration purpose. Fig. |7Ja) shows an anti-correlated case, in which both silicon/oxide interfaces shift 
outward to make a thicker body. Fig. [TJb) has correlated variations, making the channel width constant. Fig. [TJc) 
also has an anti-correlated case but both interfaces shift inward, making a thinner body in the middle of channel. 
For each case, oxide volume change compensates for silicon volume change and vice versa, so that all of them 
appear to be the case of silicon/oxide interface displacements. The change of drain current by applying 1% of 
geometric variation from numerical experiments (markers) and from GV (lines) are compared in Fig. [8ja). We can 
see that GV model’s predictions agree well with the results of numerical experiments. 

A/d-Vg results from GV model in Fig. [8] should be combined with the Id~Vg of the nominal device to recover 
the Id-Vq curves of all variational devices. Since all variations are from the same nominal device, it is convenient to 
evaluate the performance of variational devices on the AI d -Vq curves directly. Among the three cases, anti-correlated 
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Fig. 7. Illustration of geometric variation types, (a), (b), and (c) are interface displacements; (a’) (b’), and (c’) are body deformation, (a) and 
(a’) are anti-correlated outward; (b) and (b’) are correlated; (c) and (c’) are anti-correlated inward. 



Fig. 8. Comparison of drain current change from numerical experiments (markers) and GV model (lines) for (a) silicon interface displacements 
and (b) body deformation for Fig. [7] All values are positive except cross symbols (x) whose absolute values are used instead in the log-scale 
plot. Vc=0.8 V. 


outward shape (black circle) corresponding to Fig. [TJa ) outperforms the other two regarding both the sub-threshold 
slope and the on-current because thinner oxide thickness provides a better gate control. On the contrary, thicker 
oxide on both gates in the anti-correlated inward shape (red diamond) corresponding to Fig. [TJc) performs the worst 
(cross symbols stand for negative values). 

2) Body Deformation: Three body deformation cases (a’), (b’), and (c’) in Fig.[7]are derived from (a), (b), and (c) 
by applying the same silicon variational shapes, except keeping the oxide thickness unchanged. Fig. [8jb) compares 
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Fig. 9. Three dimensional view of a half structure nMOS FinFET on a SiC >2 substrate. Source/drain doping is of 10 20 cm 3 and the channel 
is intrinsic. 


the results of GV model (lines) with numerical experiments (markers), which also show good agreement. It also 
reveals that the anti-correlated inward case (red diamond) for Fig. [TJc’) has the best electrostatic characteristics 
among the three: its on-current increases and its sub-threshold slope becomes steeper because of better gate control 
on the thin body. 

For both silicon/oxide interface displacement and body deformation cases, GV model has shown accurate 
predictions confirmed by numerical experiments. Other geometric variation situations may be considered involving 
other body deformations which GV model can be easily applied to. Another benefit of GV model is that it can 
obtain terminal characteristic changes for all samples during the one nominal device simulation, while numerical 
experiments require running each sample separately. 


B. Three Dimensional FinFET 

For a more realistic application, we use our GV model to study the geometric variations of an n-type FinFET as 
shown in Fig. [9] The FinFET has the channel length of 15 nm, the body thickness of 6 nm (half structure of 3 nm 
is shown), and the fin height of 15 nm. In the following simulations, we keep drain voltage at 0.8V. Gate voltage 
increases from 0 to 0.8V. 

1) Linearity: Multiple variation amplitudes from 2% to 10% with a steady step of 2% have been applied uniformly 
on the device channel length, body width, and fin height. Corresponding drain current changes are calculated using 


the GV model, which are compared with numerical experiments in Fig. 10 Generally, GV results match with 
numerical experiments and the accuracy decreases with the increase of noise source magnitude. Comparing the 


absolute values of A/d in the three cases, Fig. 10 also reveals that given the same geometrical variation ratio, the 
on-current of the FinFET in Fig. [9] is most sensitive to the height variation (Fig. [T0|c)), followed by the width 
(Fig. [TOfo)), and least sensitive to the length variation (Fig. [I0|a)). With an increase of channel length, it is shown 
that drain current decreases (negative Aid). In comparison, the expanded cross section by the increase of the body 
width or fin height leads to the increase of the drain current. 
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Fig. 10. A/d-Vq curves with 2% increment of geometric variation from 2% to 10% on (a) length, (b) width, and (c) height. Arrows point to 
the increase of variations. 

2) Effective Noise Source: As discussed in the effective noise source calculation, calculating SAT has to take into 
account not only the box coefficient and volume changes, but also the change of flux and source terms in 0. For 
the physical models we have applied, four components are considered for the geometric variations: box coefficient 
A and volume Cl (Box coef), node distance to the interface (Dist intf), gradient of box coefficient (Grad coef), 
and quantum potential boundary condition (QPot be). Fig. El shows the each noise source contribution by the 2% 
variation of length, width, and height. Only considering box coefficient and volume (Box coef), GV’s predictions 
can be quite off in the case of length and width variations. In the length variation (Fig. [Tlja)), node distance to 
interfaces does not change so that it (Dist intf) gives a zero influence on the total result. Fig.[lTJb) suggests that Box 
coef, Dist intf, and Grad coef contribute to the same direction in the width variation. Height variation is dominated 
by the box coefficient and volume (Box coef), and net effect of others only provides minor corrections. 

3) Equation Contribution: For the case of the 2% variation, a detailed examination of the contribution from 
each equation in 0-0 for the three cases is illustrated in Fig. [12] In the length variation (Fig. [l2ja)), the electron 
continuity and Poisson’s equations are two main contributing sources, and both of them decrease the drain current. 
Increasing either body width (Fig. [f^b)) or fin height (Fig. [l2jc)) enlarges the cross section, so that the electron 
continuity equation brings up the drain current. An enlarged cross section also reduces the quantum confinement, 
thus electron density-gradient potential equation contributes a small positive amount. The negative contribution of 
Poisson’s equation in the width case (Fig.[T2|b)) cancels out a significant part of the total drain current change, while 
in the height variation (Fig. [l2|c)) the effect of Poisson’s equation is small, only comparable to the density-gradient 
equation’s contribution, leaving the net drain current change larger in the height variation case. 
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Fig. 11. Contributions of effective noise source components on drain current change for 2% geometric variations of (a) length, (b) width, and 
(c) height. 



Fig. 12. Individual equation and their total contribution on drain current change for 2% geometric variations of (a) length, (b) width, and (c) 
height. For nMOS, hole contributions from the continuity and density-gradient quantum potential equations are negligible for all cases so they 
are not shown. 


4) Sensitivity Map: The last benefit of applying GV model is that geometrically sensitive areas which mostly 
influence terminal characteristics can be easily identified. At both drain and gate voltages equal to 0.8V, for the 
three variational cases we can plot sensitivity maps. Fig. [13] shows a plane along transport direction close to the 
SiC >2 substrate. The plots show the node-wise absolute value of local noise of the drain current in the log scale for 
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Fig. 13. Sensitivity maps of equation contributions of drain current local noise for 2% geometric variations of channel length. 


a uniform 2% channel length variation. We can see that the electron density gradient equation gives contribution 


(Fig. [L3[c)) of two orders magnitude smaller than the Poisson’s and electron continuity equations. The largest 
drain current local spatial noise for Poisson’s equation comes from the source/channel boundary where we have 
applied geometric variation and also the channel/oxide interface close to the source side (Fig. [l3ja)). For the 
electron continuity equation, both source/channel and drain/channel boundaries are important although source side 


is contributing more (Fig. 13 T>)). Similarly for the body width variation, sensitivity maps are drawn in Fig. 14 


Source/channel overlap region is important for the Poisson’s equation (Fig. [I4ja)). A small part of drain region 
which is adjacent to the channel shows high noise from the electron continuity equation (Fig. [l4jb)). For height 
variation, Fig. |T5jb) indicates that the leading contribution from the electron continuity equation comes from the 
drain/channel interface. Poisson’s equation’s contribution mainly comes from the silicon/oxide corners on the source 
side (Fig. [I5ja)), which is similar to the case of the electron density gradient equation (Fig. [I5{c)). The quality of 
source/channel interface also has a big impact on the noise contribution from the electron density gradient equation 
(Fig.[l5jc)). One should note that because nonlocal effects (distance to interface, mesh gradient, etc.) are considered 
in the local noise source calculation, the plotted local contribution in Figs. [T3p~5| appears to be somewhat noisy 
especially for the most affected electron continuity equation case. In all, Poisson’s and electron density gradient 
equations are more sensitive to the source/channel interface, while electron continuity equation depends on the 
drain/channel interface. 


V. Conclusions 

To accurately and efficiently address the increasing importance of geometric variations in nanoscale electronic 
devices, we have proposed the GV model. As a comparison, we have explored physical limitations of the existing 
RGF model regarding interface displacement and body deformation. The fundamental difference between the 
proposed GV model and RGF model is on the noise source origins. GV considers the geometric variation directly 
by mesh deformation, and computes the corresponding noise sources consistently with the mesh discretization and 
the employed physical models. On the contrary, RGF model relies only on the interface quantities to compute the 
effective noise source. As a result, the RGF model cannot capture the important nonlocal geometric dependences. 
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Fig. 14. Sensitivity maps of equation contributions of drain current local noise for 2% geometric variations of body width. 
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Fig. 15. Sensitivity maps of equation contributions of drain current local noise for 2% geometric variations of fin height. 


GV, on the other hand, relies on the Green’s functions of coupled equations which are intrinsically globally defined, 
and thus applicable to wider variety of geometric variations. 

Regardless of the interface displacement or body deformation, by taking into account all effective noise source 
components: box coefficient and volume, distance to the interface, gradient of box coefficient, and proper boundary 
condition, GV model can accurately predict the changes of terminal characteristics, which have been extensively 
validated by comparing with the results from numerical experiments of full device simulations. 

GV model also offers great details for internal quantities: contributions of effective noise source components 
can be compared quantitatively, equation-wise noise contributions can be also mutually compared, and node-wise 
geometric sensitivity maps can be used for identifying sensitive regions. 

GV model can be a useful tool for geometric variability analysis. 
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